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Confocal microscopy of colloids combined with digital image processing has become a powerful tool 
in soft matter physics and materials science. Together, these techniques enable locating and tracking 
of more than half a million individual colloidal particles at once. However, despite improvements in 
locating algorithms that improve position accuracy, it remains challenging to locate all particles in 
a densely-packed, three dimensional colloid without erroneously identifying the same particle more 
than once. We present a simple, iterative algorithm that mitigates both the “missed particle” and 
“double counting” problems while simultaneously reducing sensitivity to the specific choice of input 
parameters. It is also useful for analyzing images with spatially-varying brightness in which a single 
set of input parameters is not appropriate for all particles. The algorithm is easy to implement and 
compatible with existing particle locating software. 


Colloids are materials comprised of nanometer- to 
micrometer-scale solid particles suspended in a fluid. 
Colloidal particles can easily be made to assemble into 
a variety of structures and phases, including gels, crys¬ 
tals, and glasses m, which can then be manipulated 
locally or as macroscopic materials. Three dimensional 
confocal microscopy of colloids combined with fast com¬ 
puters and precise image processing makes it possible to 
track individual particle locations over time in a bulk 
colloid mg. These techniques enable precise measure¬ 
ment of the structure and dynamics of materials over the 
entire range of length scales from the constituent parti¬ 
cles up to the bulk. They have been used for investiga¬ 
tions into colloid structure formation EHin], gel structure 
and dynamics EHni, dynamical heterogeneities and lo¬ 
cal deformation mechanisms in glasses nnsHn], aging of 
glasses [18] , and defects and elasticity in crystals [IQIIIS- 
[ 22 |. 

However, the information to be gained from these ex¬ 
periments is only as accurate as the particle locations 
themselves, since they are the foundation upon which all 
further analyses are based. For this reason, it is essen¬ 
tial that the particles be located as precisely, accurately, 
and completely as possible. There exist several general 
approaches to object detection in images [23j [24], but 
these methods are usually not accurate enough for col¬ 
loidal particle location [25]. Standard algorithms do ex¬ 
ist for the specific problem of precisely locating spherical 
objects in 3D images, and which address the particu¬ 
lar problem of finding particles in colloidal suspensions 
[26] . Several recent studies have focused on algorithms to 
improve the accuracy and precision of particle locating 
[27H3Q] . However, significant challenges remain, particu¬ 
larly in analyzing 3D images of densely-packed particles 
or images with spatially-varying brightness. While indi¬ 
vidual particles may be located very accurately, it is very 


common for some particles either to be missed entirely 
or identified more than once (“double-counted”). A user 
can avoid double-counted particles either by enforcing a 
strict minimum separation between particle locations or 
by adjusting input parameters to lower the sensitivity of 
the locating software, but at the cost of missing more 
particles. On the other hand, increasing the software 
sensitivity also increases the number of particles found 
multiple times; although setting a separation threshold 
works fairly well to limit this double-counting, it intro¬ 
duces an additional parameter that may strongly affect 
final results and there is still no guarantee that all parti¬ 
cles will be found. The result is a tradeoff between these 
two types of errors that is very sensitive to the user’s 
precise choice of input parameters. Attempts to choose 
input parameters that mitigate both problems simulta¬ 
neously inevitably produce some of each type of error. 

To solve this problem, we developed an iterative parti¬ 
cle locating algorithm that is able to find all of the par¬ 
ticles in a 3D image of a colloid with few or no double- 
counted particles. We first run standard particle locating 
software on the original image data, using input param¬ 
eters that deliberately err on the side of missing parti¬ 
cles in order to avoid double-counted particles. Next, we 
erase from the original image those particles that have 
already been found. This creates a new “raw” residual 
image that is empty except for those particles that have 
not yet been found. We then apply the particle loca¬ 
tion software to the residual raw image, and iterate this 
procedure until all particles have been found. 

The steps of the iterative algorithm are described in de¬ 
tail below, illustrated for an example image stack in Fig¬ 
ure For a typical 3D image containing 50,000 particles 
of equal brightness, we locate about 95% of the particles 
on the initial pass through the data, and the iterative 
locating is complete after four to six iterations. At the 
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end, we estimate that <0.2% of the particles are double- 
counted, and these isolated pairs can easily be identified 
and consolidated. For a given sample, algorithm perfor¬ 
mance will depend on the resolution and quality of the 
original image data; for example, images with spatially- 
varying brightness, such as the example shown in Figure 
often require more iterations than evenly-illuminated 
samples. 
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FIG. 1. Iterative particle locating applied to a 3D confocal 
image of a dense, disordered colloid with variation in indi¬ 
vidual particle brightness and uneven illumination across the 
sample ED. (a) An x-y cross section through the raw data. 

(b) The same cross section through the residual raw data af¬ 
ter particles found in the first locating pass have been deleted. 

(c) The results of iterative particle locating. Particles whose 
centers are within about ±1 particle radius of this cross sec¬ 
tion are marked according to the locating iteration in which 
they were found: 1st pass, (o); 2nd pass, (A); 3rd pass, (o); 
4th pass, (□); 5th pass, (y)- (d) The results for the same 
3D data set of standard non-iterative particle locating with 
optimized parameters. 


The software we use to implement iterative particle lo¬ 
cating is included as MATLAB code in the Supplement 
ED- We also provide three example data sets: Example 
1, a colloidal glass; Example 2, a colloidal crystal; and 
Example 3, a dense, disordered colloid with spatially- 
varying image brightness (used for the example shown in 
EigureQ. Each example data set includes a small 3D 
confocal image stack, a text file of particle locations, and 
the input parameters used to locate the particles. We in¬ 
terface the iterative algorithm with the publicly-available 
MATLAB particle locating software of Ref. m- How¬ 
ever, the algorithm we present is broadly compatible with 
any approach to particle locating, and the example code 
provided should be easily adaptable to interface with 
other locating software. 

We start with a 3D confocal image of a colloid in which 
the particles appear bright against a dark background. 
An x-y cross section through such an image is shown in 
Eigure[^a). We bandpass filter the raw image and run a 
first pass of particle locating using the standard filtering 
and feature-finding software of Ref. [29] . In choosing in¬ 
put parameters for the standard software, we deliberately 
choose parameters that may miss some particles rather 
than double-counting any particle. Any missed particles 
will be identified in subsequent iterations. 

Next, we generate a 3D residual raw image by deleting 
from the original image those particles that have already 
been found. Spherical particles generally appear as ellip¬ 
soids in 3D images, so we erase each found particle by 
replacing the original image data surrounding its coordi¬ 
nates with an appropriately-sized ellipsoid of zero-valued 
pixels. The x, y, and z particle dimensions in pixels are 
the only input parameters required by the iterative lo¬ 
cating algorithm beyond standard particle locating. 

The resulting residual raw image contains only those 
particles that have not yet been found, as shown for the 
example cross section in Eigurej^b). In this case, a num¬ 
ber of particles were missed during the first pass. Parti¬ 
cles that were missed are usually isolated in the residual 
raw image, as on the left side of the example, making 
them much easier to locate on subsequent iterations than 
when they were surrounded by other bright objects. In 
cases where spatially-varying image brightness causes the 
initial locating pass to miss entire regions of particles, as 
in the upper-right corner of the example, the iterative 
locating procedure will locate all of the particles over the 
course of several iterations. 

We then bandpass filter the residual raw image and run 
the particle locating software again. We use all the same 
input parameters for filtering and locating as in the first 
pass, with the added constraint that a bright region must 
have an integrated intensity greater than some minimum 
threshold in order to be considered a real particle rather 
than noise. This is a standard feature of many particle 
locating software packages, including the one we use [29] . 
This cutoff prevents any small bright regions in the resid- 
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ual raw image from being falsely identified as additional 
particles, such as the edges of particles that may have 
been incompletely deleted from the original raw image. 
As these regions are significantly smaller than the real 
particles, there is no ambiguity in distinguishing them. 
Particles cut off by the edges of the image may or may 
not pass the integrated intensity threshold. Although 
they are real particles, their center coordinates will not 
be accurate if they are not fully contained within the 
image. Fortunately, it is straightforward to exclude par¬ 
ticles close to the image edges from subsequent analyses. 

We continue to delete particles from the images as they 
are located, and iterate this entire process until no new 
particles are found. For a typical image with uniform 
brightness, the second iteration finds nearly all of the 
particles missed during the first pass, and the locating is 
usually complete in 4-6 iterations. Samples that were un¬ 
evenly illuminated by the microscope or that have a large 
variation in the brightness of individual particles may re¬ 
quire more iterations. These are also data sets for which 
a single set of locating parameters would usually not suc¬ 
cessfully identify all particles, so the iterative approach 
is particularly helpful here. The 3D data set shown in 
the Figure required 8 iterations to complete the locating. 

Figure [^c) shows all of the particles located near this 
cross-section at the end of the iterative locating process. 
Each found particle is marked according to the iteration 
in which it was located. For comparison, we also show 
the results of traditional single-pass particle locating in 
Figure [^d). In this case, we optimized the parameters 
to locate as many particles as possible in a single pass 
with minimal double-counting, including setting a mini¬ 
mum separation distance between neighboring particles. 
However, the locating results are only marginally more 
complete than the first pass of the iterative particle lo¬ 
cating, and many particles in this example were missed 
entirely. 

In summary, we have created and implemented an 
iterative algorithm to improve the completeness parti¬ 
cle locating of individual colloidal particles in an image. 
The algorithm is useful, simple, straightforward to im¬ 
plement, and easily integrated with existing particle lo¬ 
cating software. The algorithm requires only three new 
input parameters from the user: the x, y, and z dimen¬ 
sions in pixels of an individual particle in the original 
image. Although we implement it for locating monodis- 
perse, spherical particles, this algorithm could easily be 
extended to colloids comprised of polydisperse or non- 
spherical particles. In the latter case, the orientation 
of the already-found particles would also be required to 
erase them accurately from the raw images. 
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